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The objective of this research is to study small transition metal clusters 
at a high level of approximation i.e. including all the valence electrons in the 
calculation and also including extensive electron correlation. Perhaps the 
most useful end result of these studies is the qualitative information about 
the electronic structure of these small metal clusters. This information 
includes the nature of the bonding e.g. whether the 4s or 3d electrons 
are involved in the bond orbitals and also can provide information about 
reactivity. The electronic structure studies of the small clusters are directly 
applicable to problems in catalysis. From comparison of dimers, trimers 
and possibly higher clusters it is possible to extrapolate the information 
obtained to provide insights into the electronic structure of bulk transition 
metals and their interaction with other atoms and molecules at both surface 
and interior locations. 

Because of the numerous practical aplications of transition metals, 
small metal clusters are currently of considerable experimental interest, and 
some information is currently becomming available both from matrix ESR 
studies and from gas phase spectroscopy. Currently the ESR experiments 
provide information on the neutral molecules (e. g. spin state and orbital 
character of the radical orbital as well as limited geometric information 
such as whether the atoms are equivalent or inequivalent). The gas phase 
spectroscopic work in principle can provide more detailed information, but 
Tor Primers and larger clusters, only C.U3 bas been isuccesfully studied. The 
experimental difficulties arise because of the high density of electronic states 
for these molecules and because for the neutral molecules a mixture of 
clusters (i. e. monomer, dimer, and trimer etc.) are present in the molecular 
beam. Future experimental studies therefore will be made on the positive 
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ions of the metal trimers, since these molecules may be mass selected. Thus, 
the ESR experiments are being used to study the neutrals, and gas phase 
spectroscopy is expected to provide information on the positive ions; while 
theory can provide a link between these experiments, since it can study 
either neutrals or positive ions. 

The studies in the past year have concentrated on the electronic struc- 
ture of the transition metal trimers. Here detailed studies have been made 
of Sc 3 and Cu 3 which are the elements on the ends of the first transition 
row. Studies are also being carried out for the positive ions of the transition 
metal trimers to help the experimentalists in their gas phase spectroscopic 
studies. Calculations have also been started for Ti^ - and V^ - and these will 
be continued in the current grant year. Comparison of these calculations 
to those for ScJ~ provides very valuable insights into the bonding trends 
in these molecules, just as similiar studies for the transition metal dimers 
provided a simple picture of the electronic structure of these molecules. As 
a start in understanding the reactivity of these molecules, the interaction 
of Sc 3 with a H atom will be studied. Initially, very simple questions such 
.ns whether the H atom will bond to the 3d or 4s derived levels of Sc 3 will 
be considered. 

The studies on reactivity are relevant to experimental studies by Smalley 
et. al. on the reactivity of H 2 with Nb clusters and studies at Argonne 
National Laboratory on the reactivity of Fe clusters. 

The main topic of this report is the electronic -structure studies for Sc 3 
and Cu 3 . Preprints of this work .are included in the appendix and only a 
few of the major results are presented here. 

For the Cus molecule all-electron studies have been carried out for 


2 



* 


06 


the 2 E state arising from three 4s 1 3d 10 Cu atoms. For equilateral triangle 
geometries(D 3 / l ) the lowest state is found to be 2 E' arising from 4sa , 2 4se' 1 . 
The 2 E' state exhibits strong Jahn- Teller distortion, leading to 2 Ai (acute 
angle) and 2 B 2 (obtuse angle) minima in C 2 V symmetry. Here the 2 B 2 
minimum is a true minimum on the surface while the 2 Ai minimum is 
a saddle point connecting adjacent 2 B 2 minima. This strong Jahn- Teller 
distortion is consistent with the observed ground state vibrational levels in 
the gas phase spectroscopic studies of Morse, Hopkins, Langridge-Smith, 
and Smalley. The 2 B 2 minimum is also consistent with the observed ESR 
spectrum of CU 3 in a matrix, which has been interpreted as an obtuse angle 
structure with most of the spin density on the end Cu atoms. A linear 2 £+ 
state is found to be 0.26 eV higher. 

Calculated vertical excitation energies for C113 indicate that the upper 
state in the transition studied by Morse et. al. should be 2 A' l not 2 E". 
An alternative assignment based on a 2 A , 1 upper state is shown to account 
for the pattern of allowed hot bands and this revised assignment has been 
accepted by the experimentalists. 

Calculations have also being carried out for Sc 3 molecule. This system 
was selected for study because of the ESR studies of Weltner et.al. which 
establish that the ground state of Sc 3 is a doublet state with equivalent 
Sc atoms which implies pairing of the 3d electrons into bond pairs. The 
.calculations .are in Agreement with experiment in that the doublet states are 
lower than the high spin states and the wavefunctions show very definite 
indications of 3d bonding. Three doublet states bave been studied. The 
dominant configurations for these states are: 

2 Aj = 4saa , 1 2 4 s< 7 e /4 Sdcra'/ 3 d 7 ra 2 /2 
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2 A 2 = 4saa , 1 2 4scre' 4 Sdtxa' 2 Sdno!^ 1 

2 E' = 4soa , 1 2 4s<je /3 3d(7a' 2 3d7ra2 2 

Based on a combination of the calculations and the ESR experiments 
it was concluded that the 2 A^ state is the ground state of Sc 3 . Calculations 
were also carried out for the Sc^ - molecule in order to assist the experimen- 
talists in doing gas phase studies of this species. 

The results of this work will be presented at several scientific meetings. 
Among these are a i) NATO workshop on ”Ab-Initio Quantum Chemistry- 
The Challenge of Transition Metals”, to be held in Strasbourg, France, Sept. 
1985 (invited 30 minute talk) ii) 5th International Congress on Quantum 
Chemistry, to be held in Montreal, Canada (poster presentation), and iii) 
an ACS Symposium on ” The Theory of Metal Metal Bonding” to be held 
in New York city, April 1986 (Invited Talk). 
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Abstract 

CASSCF/CCI calculations are presented for the low-lying states of SC3 
and Sc^ - , and SCF/CI calculations are presented for the 1 A , 1 state of Ca 3 
arising from three ground state (4s 2 ) Ca atoms, all for equilateral triangle 
geometries. The calculations use effective core potentials, developed by Hay 
and Wadt, which replace the Ne core but include the 3s and 3p core levels 
along with the valence electrons in the calculations. 

The bonding in Ca 3 arises by 4s-*4p promotion and leads to a well 
depth of about 0.5 eV for R(Ca-Ca) = 7.5 ao. For Sc 3 the 4s bonding 
is simiiiar to that in Ca3, but the 3d electrons are also strongly bonding 
leading to a 2 A2 ground state with a well depth of about 1.0 eV and R(Sc-Sc) 
— 5.75 ao- The good 3d bonding orbitals (bonding between all three atoms) 
are 3da^ derived from atomic 3d7r" and 3da' x derived from atomic 3dcr, while 
Sdn^ atomic orbitals lead to 3de' orbitals which are bonding between pairs 
of atoms, and the 3d5' and 3d$" derived levels are non-bonding. (Here the 
atomic symmetry is given with respect to an axis connecting the atom to 
the center of the molecule.) Based on the Sc 3 calculations and preliminary 
calculations on Ti^", the bonding in V^ - and Cr^ - is also discussed. 


a Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 



I. Introduction 


Transition metals and transition metal(TM) compounds are currently 
of considerable interest because of their relevance to catalysis and to materials 
science problems such as hydrogen embrittlement and crack propagation in 
metals. In spite of such interest, progress in understanding the chemistry of 
TM molecules has been slow due to both experimental and theoretical 
difficulties. For a review of the current state of both theory and experiment 
for TM molecules the reader is referred to the review article by Weltner and 
Van Zee[l]. 

For the TM trimers most of the current experimental information 
comes from matrix ESR studies. Here the assumption is made that the 
electronic state observed in the matrix is the ground state of the molecule. 
These studies directly give the electronic spin state of the molecule and 
analysis of the hyperfine structure gives some information on the electronic 
distribution of the open shell orbital(s). In some cases limited geometric 
information can also be derived e.g. for Sc 3 all the Sc atoms are found to 
be equivalent thus implying an equilateral triangle geometry[2]. To date 
gas phase spectra have been obtained only for the copper trimer[3]. These 
spectra permit only a vibrational analysis and thus do not give detailed 
geometric information, but the vibrational analysis does indicate a Jahn- 
Teller distorted ground state and provides vibrational frequencies for an 
excited state at an excitation energy of about 2.3 eV. 

Because of the large number of low-lying states expected for the TM 
trimers and because of the presence of other multimers, the resonant two- 
photon ionization experiments, which were successful for CU 3 , can not 
be used in general to study the remaining TM trimers, due to spectral 
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complexity. It is, however, possible to use mass spectroscopy to pre-select 
and focus the positive ions. Thus, it is likely that the gas phase work 
will yield information about the positive ions, while ESR work will focus 
on the neutrals. Theory can provide a link between these two classes of 
experiments, since it can be applied to both the neutrals and positive ions. 

In this paper we discuss complete active space SCF[4]/ externally 
contracted[5] configuration interaction (CASSCF/CCI) calculations for the 
low-lying states of SC 3 and Sc^\ We also compare the bonding in the 
Ca 3 , SC 3 , and Cus[ 6 ] molecules. From this comparison we are able to 
predict general trends for the TM trimers, much as we observed for the 
TM dimers [7]. 

Section II discusses the qualitative features of the bonding in the TM 
trimers including predicted trends for other TM trimers. Section IH dis- 
cusses the basis sets and other technical details of the calculations. Section 
IV discusses the calculated results for Sc 3 and Sc^ and compares the results 
to experiment. Finally, Section V presents the conclusions from this work. 
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II. Qualitative Discussion. 

As with the TM dimers[7] the nature of the bonding in the TM trimers 
changes with position in the transition row. The two principle factors 
controlling the electronic structure are i)the relative (n-fl)s and nd orbital 
sizes which control the contribution of (n-j-l)s and nd bonding and ii) 
the (n4“l)s 2 nd m -*■ (n-f-l^nd” 1 "*" 1 excitation energy of the atom which 
controls the atomic configurations which are available for bonding. As 
discussed in ref. 7, the ratio <T( n +i)6> / <Tnd> increases monotonically 
from 2.03 for Sc to 3.36 for Cu. Thus, we expect that nd bonding will be 
more favorable for the Sc end of the row than for the Cu end of the row. 
For Sc the 4s 2 3d 1 state of the atom is lowest with the 4s 1 3d 2 state about 
1.4 eV higher. This excitation energy decreases monotonically from Sc to 
Cr which has a 4s 1 3d 5 ground state. For Mn, 4s 2 3d 5 is again lowest due to 
the large exchange interaction associated with a half filled 3d shell, but as 
for the left half of the row, the excitation energy decreases monotonically 
from Mn to Cu, which has a 4s 1 3d 10 ground state. Thus, we expect that 
the ground states for Cr 3 and Cu 3 will arise from the 4s 1 3d m+1 states, 
while for the Sc 3 molecule the ground state configuration will involve a 
compromise between possible better bonding arising from 4s-+3d or 4s-+4p 
promotion and the excitation energy involved. Moving from Sc 3 toward Cr 3 
one expects more 4s-»-3d promotion as the (n+l)s 2 nd m -*■ (n-f-l)s 1 nd f ”+ 1 
excitation energy decreases. 

For Cus the ground state of the atom is 4s 1 3d 10 with the closed 3d 
shell contracted with respect to the 4s. Thus, the bonding in Cus is found 
to involve mainly the 4s electrons[6], much like an alkali trimer. However, 
the 3d electrons are not inert since correlating them results in a decrease 
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in R Cu-Cu of about 0.30 ao- The Mulliken populations for CU 3 from 
a Cl calculation at the optimal geometry for Dqh symmetry (the ground 
state geometry here is Jahn- Teller distorted but near equilateral triangle) 
are 4s=2.72, 4p=0.47, and 3d=29.80 in agreement with the qualitative 
picture above. 

For Ca 3 the lowest atomic limit is for three ground state atoms ( 
*S or 4s 2 ). SCF/CI calculations have been carried out for this limit for 
equilateral triangle geometries and lead to the potential curves given in 
Fig. 1. Here the SCF curve is repulsive, but the Cl curve shows a weak 
bond (about 0.5 eV) at R Ca-Ca = 7.5 ao- At this distance, the populations 
from the Cl calculation are 4s=5.24, 4p=0.61, and 3d=0.15. At R Ca-Ca 
= 6.5 ao, the populations are 4s=4.58, 4p=1.08, and 3d=0.33. These 
results indicate a significant amount of 4s -+ 4p promotion as well as some 
4s-*- 3d promotion for shorter R Ca-Ca. For Ca 3 this formation of sp 
hybrids apparently does not lead to strong bonding; however, the large 4p 
populations at smaller R indicate that the predominate bonding mechanism 
here is 4s-*- 4p promotion, as compared to the C &2 molecule which does not 
show the large 4p population and where a large part of the bonding is due 
to Van der Waals terms arising out of the 4s-*- 4p near degeneracy effect[ 8 ]. 
The larger importance of 4s -*■ 4p promotion for the trimer as compared 
to the dimer probably results from the formation of three bonds for the 
trimer which compensates for the 4s -*■ 4p promotion energy, while the 
dimer with only one bond can not compensate for the promotion energy. 
Similar results have been observed by Bauschlicher, Bagus, and Cox[9] for 
Be 4 , Mg 4 , and Ca 4 . 

For SC 3 we find the low-lying states arise from two atomic limits i) 
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three ground state 4s 2 3d 1 atoms and ii) two ground state 4s 2 3d J atoms 
and one excited state 4s 1 3d 2 atom. This is to be expected based upon the 
results for Sc2[10],[ll], 'where it was shown that the ground state arises 
from one atom in the 4s 2 3d 1 state and one atom in the 4s 1 3d 2 state. For 
both of these atomic limits the bonding involves 4s to 4p promotion leading 
“to sp hybrid bonding from the 4s shell, in addition li) also involves 4s-+3d 
promotion leading to a 3d population of about 4.0. For Sc 3 the 3d orbitals 
also contribute to the bonding. Because of the smaller size of the 3d orbitals 
compared to the 4s orbitals “the formation of 3d bonds results in a much 
shorter Sc-Sc distance in Sc 3 as compared to the CarCa distance in Ca 3 . 
This is shown in Figure 2 where one sees that the 2 E ; state of Sc 3 has a 
minimum, for D 3 * geometries, for R Sc-Sc= 5.5 ao, while the Ca 3 potential 
curve is strongly repulsive at this distance. This result implies a strong 
bonding effect from the 3d shell; although, the size mismatch of the 4s 
and 3d orbitals results in repulsive 4s-4s interactions at the short R Sc-Sc 
distances needed for formation of good 3d-3d bonds, and results in only 
“weak molecular bonding. 

The two principle 3d bonding orbitals for Sc 3 are a^ “which arises from 
the 3d<r atomic orbitals and a^ "which arises from the SdTr" atomic orbitals. 
The configurations of the low-lying states here are: 

^A' = 4sa^ 2 4se' 4 3da' 1 3da" 2 

s 4sa'j 2 4se' 4 3da^ 2 3da^' 1 
— 4sa 4se ,s 3da^ 2 3da^ 2 

The populations, for CASSCF wavefunctions at equilateral triangle geometries 
with R Sc-Sc= 5.5 ao, are: 4s=4.55, 4p=1.17, 3d=3.24, and 4f=0.04 for 
the 2 A^ and 4s=3.89, 4p=1.24, 3d=3.83, and 4f=0.04 for the 2 E'. Thus, 
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the wavefunctions contain a considerable admixture of 4p character in 
agreement with the discusssion above. 

The 2 Aj and 2 A% states, which arise from three ground state atoms 
are weakly repulsive at the CASSCF level, as for Ca 3 , but are bound at 
the Cl level. The 2 E' state, which arises from two ground state atoms and 
one excited state atom(4s 1 3d 2 ), is bound at the CASSCF level only with 
respect to its asymptote but is also bound with respect to three ground state 
atoms at the Cl level. The 2 E' state is calculated to be 0.1 eV above the 
2 A 2 state, in agreement with the ESR results[2] which indicate a non Jahn- 
Teller distorted ground state, consistent with the calculated 2 A 2 ground 
state of SC 3 . 

For Sc+, the promotion energy 4s 1 3d 1 -*• 3d 2 is much less than 4s 2 3d x 
-*• 4s 1 3d 2 for Sc atom [12], thus, for Sc^ one expects a higher 3d population 
to be favored, and the ground state of the ion is expected to be derived 
from the 2 E' state of the neutral. By analogy to CU 3 , where the 4se' orbital 
is only weakly bound, the lowest state of Sc^ arises by removal of one 
electron from the 4se' orbital. This leads to a 3 A^ state with the electronic 
configuration: 

8 Ag = 4sa , 1 2 4se ' 2 3da , 1 2 Sda ^ 2 

This configuration also leads to higher-lying 1 E' and 1 A / 1 states. Starting 
with the 2 Aj, and 2 A^ states, removal of one electron from the 4se' orbital 
leads to S E' and 3 E ,/ states with the configurations: 

3 E' = 4sa ' 2 4se /3 3ddk [ 1 3 da 2 2 
8 E" = 4sai 2 4se ' 3 3daj 2 3da^ ,! 

These states are at 0.58 eV and 0.66 eV respectively above the ground state 




We now consider what would be likely configurations for other TM 
trimers of the first transition row. As Z increases, there will not be enough 
good 3d bonding orbitals available to accomodate all the 3d electrons. This 
situation is analagous to the TM dimers, where the 3dcr ff and 3d7r u are 
good bonding orbitals and the 3d<5 tf are weak bonding orbitals. Thus, 
the TM dimers are able to accomodate six electrons in good bond or- 
bitals. For Ti 2 the 3d orbital occupancy is 3d<7 2 3d7r*. For V 2 the ex- 
tra two electrons go into the weakly bonding 3d<5 ? orbitals leading to the 
configuration 3dcr 2 3d7r*3d<5 2 . Here the 6 orbitals form one-electron bonds 
which do not contribute significantly to the bond strength. For Cr 2 the 
3d 6 g orbitals contain four electrons which is unfavorable and the binding 
energy in Cr 2 is less than in Y 2 . This is understandable since for weakly 
overlapping orbitals, one-electron bonds whose bonding terms vary with 
distance like the overlap(S) are more favorable than two electron bonds 
"whose bonding terms vary with distance like S 2 [10]. 

hJext we consider the symmetries of the molecular orbitals derived 
from the five possible atomic 3d orbitals. Here we consider 3d functions on 
oach atom where the sigma direction Is along the axis connecting the given 
atom with the center of the :molecule. These 3d functions are classified 
with respect to the above axis as 3d<7, Sd^, 3d7r", 3d5', and 3d<5 ,/ , where 
the primes denote symmetry with respect to the molecular plane. In the 
anolecular symmetry 3dcrand 3d<5' transform as a^ and o^Sd^ transforms 
as ag and -e^d^' transforms as a^' and e'^ and 3d$ ;/ transforms as a" and 
e". In order to -show the symmetry properties of the 3d derived orbitals, 
plots of molecular orbitals for the Ti 3 molecule arising from 3d<r, Sd^, 
3d7r", and 3d 5' are given in Fig. 3. As for the TM dimers, we find that 
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the 3dcr, 3(37^, and 3d it" atomic orbitals have moderate overlaps in the 
molecule and form the bonding orbitals, while the 3d6' and 3d5" orbitals 
have small overlaps and are essentially non-bonding. Here we see that the 
best bond orbitals are Sdaa,^ and Sd^'a^, because these combinations of 
the atomic orbitals have no molecular nodal planes between pairs of atoms. 
These same atomic orbitals also lead to 3dae' and Sd^'e" orbitals which 
contain one molecular nodal plane, and are thus not favorable as bonding 
orbitals. However, these orbitals are important correlating orbitals, since 
they introduce molecular nodal planes between pairs of atoms. The 
orbital is anti-bonding with modal planes between all pairs of atoms, while 
the 3d7rV orbitals contain one molecular nodal plane and are thus less 
favorable than Sdcra'j or 3d7r", but are still bonding between pairs of atoms. 
Finally, the 3dtf orbitals remain very atomic like due to the small overlaps 
of these atomic orbitals. Thus, these orbitals are essentially non-bonding. 

This ordering of the 3d molecular orbitals Is consistent with the oc- 
cupations of the low-lying -states of SC3 where the 3daa' and 3d rf' a^ orbi- 
tals are occupied. For the Ti atom, the 4s 1 3d 3 state is at 0.81 eV compared 
to 1.44 eV for the corresponding excitation in Sc. Thus, for “Ti3 one expects 
-4s-*3d promotion to be more important than for SC3. Preliminary calcula- 
tions for Ti^" Indicate two important .groups of states. The first group of 
states arises from two Ti atoms in the ground state and -one Ti+ in the 
-ds^d 2 state which .gives a 3d population of about six. These states are 
high spin (spin=8) and have an Tt(Ti-Ti) of about 7.5 ao for equilateral 
triangle geometries. The second .group of states has a 3d population of 
about seven. The wavefunctions for these states are characteristic of 3d 
bonding, and these states are stable at much shorter bondlengths R(Ti-Ti) 

t 
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= 5.0 ao- One representative configuration from this group of states is: 
1 6 Ai= 4saj 2 4se ' 2 3daa.\ 2 Sd^'a ^ 1 Sdfl^a ^ 2 3d<5V 2 

This configuration is similar to the 3 A4 state of Sc^" with the extra 3d 
electrons in 3d<5' derived levels. Since the 3d<5' levels have very small 
overlaps, they are non-bonding. An alternative configuration is: 

2 6 Aj= 4sai 2 4se ' 2 3d(jai 2 Sd^'a ^ 1 Sd^'a ^' 2 3 d 7 rV 2 

One might expect this configuration to be lower, since the 3d?rV levels are 
bonding between pairs of atoms. However, this configuration involves atomic 
configurations with two electrons in the 3d?r orbitals which is not favorable 
for the Ti atom and this intraatomic coupling effect leads to comparable 
energies for the two configurations. 

For V 3 and O 3 one expects the lowest state to arise from the three 
atoms in the 4s 1 Sd ” - *" 1 state. The low-lying configurations here are ex- 
pected to arise from the l^Aj configuration of Ti^" by first adding electrons 
to the bonding Sd^e' levels, and then adding electrons to the non-bonding 
3d5" levels. This model would result in increased bonding for V^" as com- 
pared to Ti^j", but no increase in bonding in going from to Cr^\ 
This trend is consistent with V and Cr having the largest bulk cohesive 
energies[13] for the first transition row. 

M 113 with a 4s 2 3d 5 atomic ground state is expected to be like Ca 3 with 
the bonding arising from 4s-» 4p promotion. For the remaining elements of 
the first transition row, we expect mainly 4s bonding arising from mixtures 
of the 4s 2 3d n and 4s 1 3d n+1 atomic states with the 3d electrons high spin 
coupled. Note that we expect 3d bonding to become less favorable as the 3d 
orbitals contract relative to the 4s for increasing nuclear charge within the 
first transition row. There is some experimental evidence for this viewpoint 
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from ESR studies -which sho-w that while M 112 has a singlet ground state 
but low-lying high spin states[14], the ground state of Mn 5 is high spin 
(spin=25) corresponding to all the 3d electrons high spin coupled[14]. 
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IQ. Calculational Details 

While these calculations could have been carried out at the all-electron 
level, accurate effective core potentials (ECP’s) have been developed by Hay 
and Wadt [15] which make it possible to eliminate the Ne core with very 
little loss of accuracy. Here the 3s and 3p core levels as well as the 4s, 
4p, and 3d valence levels are included in the calculation. The inclusion of 
the entire n=3 shell in the valence space is necessary since the 3p and 3d 
orbitals have similar spatial extents, e.g. removing the entire Ar core has 
lead to collapse problems for ScO[16a]. 

The basis set used in the Sc 3 calculations is the basis set reported 
by Hay and Wadt except that the 4p like functions were replaced by 
Wachters’[17] 4p functions multiplied by 1.5 to make them suitable for 
describing 4s-*- 4p correlation. The s and p basis sets were contracted 
(2111) based on the 3s and 3p atomic orbitals, respectively. Normally, one 
augments the 3d basis by the Hay diffuse 3d function[18]. However, here the 
most diffuse 3d function is 0.0812 compared to the Hay diffuse 3d of 0.0588. 
A diffuse function selected based on an even tempered criterion was tried 
but it had very little effect on the energy of either the 4s 2 3d* or 4s 1 3d 2 state 
of Sc atom and it was concluded that the valence basis set was sufficiently 
diffuse. In order to describe polarization of the 3d bonds a single set of 
4f functions was added as a two term fit of an STO(f=1.6). The Slater 
exponent here was optimized from an all electron CISD calculation for the 
4s 1 34 2 state of the atom. The final basis set is (5s5p5d2f)/[4s4p3dlf]. 

The Ca atom basis set was constructed in the same way as the Sc 
basis set and the Hay and Wadt ECP with the 3s and 3p included in the 
valence was used as for Sc. The 3d functions here were taken from the 
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work of Pettersson, Siegbahn, and Ism ail [19]. The primitive basis set here 
is 5d which is contracted (311). No 4f basis functions were added to the 
basis set, since there are no occupied 3d orbitals. The final basis set is 
(5s5p5d)/[4s4p3d). 

Several tests were carried out for the Sc basis set/ECP used here. 
Atomic SCF calculations using the ECP gave a 4s 2 3d 1 -*• 4s 1 3d 2 excita- 
tion energy of 1.01 eV compared to 1.03 eV from NHF. Also all-electron 
calculations on the 4s 1 4p 1 3d 1 state of the Sc atom showed that multiply- 
ing the 4p exponents by 1.5 raised the energy of this state by only 0.03 
eV; this result is important since we expect strong 4s-* 4p promotion in 
Sc 3 . Molecular calculations were also carried out for ScO and Sc 2 . The 
ScO results were in good agreement with the all-electron results reported 
by Pettersson, Wahlgren, and Gropen[16b]. The ECP calculations for the 
5 £“ state of Sc 2 compared well with all-electron results[20] (errors of -.04 
a 0 , -18 cm -1 , and 0.00 eV for R e , w e , and D e , respectively). 

CASSCF/CCI calculations were carried out for Sc 3 . Here the active 
space in the CASSCF included those molecular orbitals derived from 4s, 
3da, and SdTr 11 atomic orbitals. This choice is based on the discussion in 
section II. While the 4p should in principle also be included in the active 
space these functions mainly contribute to the 4s-* 4p near degeneracy effect 
which may be adequately treated in a subsequent Cl calculation. The above 
atomic orbitals lead to molecular orbitals of and e* from 4s, of and e' 
from 3d<7, and of and e" from 3d7r". This gives a CASSCF calculation of 
nine electrons in nine orbitals which is small enough that it was possible to 
do the full CASSCF calculation, i.e. no orbital occupation constraints are 
necessary. Calculations were carried out at this level in order to determine 
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the law-lying states. Subsequent calculations were then carried out with 
the number of electrons in a! and a!' orbitals constrained to be the same 
as in the dominant configuration. While this constraint results in a very 
substantial reduction in the number of configurations, it has very little effect 
on the total energy. 

The reference configurations for the CCI calculations included the 
dominant configurations involving correlation of the 3d derived levels. For 
the 2 A^ state five reference configurations were used consisting of the SCF 
plus 3daj -* e' double excitations(two configurations) and interpair terms 
between 3da'j and Sda^two configurations). For the 2 E' state seven reference 
configurations were used. These were the analagous configurations to the 
five for the 2 A" state plus the double excitations from Sda^-^de'^two 
configurations). The calculations were carried out using the MOLECULE[21]- 
SWEDEN[22] system of programs. 

IV. Discussion. 

Tables I and II give calculated CASSCF energies for SO3 and Sc^~, 
respectively, while these values are plotted in Figure 2. CCI calculations 
have also been run for the 2 A% and 2 E' states. These energies are given in 
Table HI. 

Table IV. gives SCF and SCF/CI energies for the Ca3 molecule, while 
these values are plotted in Figure 1. For the Ca3 molecule some checks were 
done on the size consistency problem. Here for D3A geometries at R Ca-Ca 
= 50.0 ao, the supermolecule Cl energy is 0.29 eV above the sum of the atom 
Cl energies, while adding Davidson’s correction[22] to the supermolecule 
Cl energy leads to an energy 0.04 eV below the sum of the atom Cl 
energies. These results indicate that Davidson’s correction works well for 
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the six electron problem resulting from three 4s 2 atoms. These results also 
indicate that the supermolecule plus Davidson’s correction energy is 'well 
approximated as the sum of the atom Cl energies and this approximation 
was used for the SC3 calculations. 

From Figure 2 one sees that the 2 E ; state has a minimum in the 
CASSCF curve near 5.5 ao. This state is derived from two ground state Sc 
atoms(4s 2 3d 1 ) and one excited state atom(4s 1 3d 2 ). The 2 E / state is bound 
with respect to this excited atomic limit but is about 0.3 eV above the 
lowest atomic limit. The 2 A£ and 2 Aj -states are slightly higher in energy 
in this Tegion and also the CASSCF curve is repulsive as ’was the case for 
the Ca3 SCF potential curve. The CCI potential curve for the 2 A^ state 
is bound by about 1.0 eV for a D3/1 geometry with R Sc-Sc= 5.75 a 0 . An 
interesting feature of the 2 A^ CCI curve is a pronounced shoulder in the 
short R region. This feature is due to a curve crossing with another state 
of the same symmetry possibly derived from the same atomic limit as the 
state. The CCI curve has a more normal shape and is calculated to 

be less than 0.1 eV above the 2 Ag state. 

The experimental ESR spectrum]2j is clearly due to "Sc3, based on the 
presence of the HFS lines and approximate intensity ratios expected for 
interaction with three equivalent Sc nuclei. Prom the magnitude of the 
nuclear HFS interaction it is concluded that the singly occupied orbital con- 
tains less than -3 % 4s character. These results are consistent with the cal- 
culated ^Ag ground state_, since this state should have an equilateral triangle 
geometry and the singly occupied orbital Is predominately "3d like. The^E' 
state should Jahn- Teller distort which would be inconsistent with equivalent 
Sc atoms. For CU3, which also has a 2 E ; ground state and exhibits strong 
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Jalm-Teller distortion[ 6 ], the lowest Jahn- Teller distorted component of 2 B 2 
symmetry is stabilized in the matrix, leading to inequivalent Cu atoms[24]. 
However, the 2 E' state of SC 3 may be only weakly Jahn-Teller distorted, 
since the 3d bonding orbitals favor equilateral triangle geometries. Thus, 
the 2 E' state could have equivalent Sc atoms if it were only weakly Jahn- 
Teller distorted and exhibits rapid pseudorotation. But even if the Sc atoms 
were equivalent, the 2 E' state would not be consistent with the ESR, since 
the singly occupied orbital is mainly of 4s-4p character. For these reasons 
~we conclude that the the 2 AJ> state is the ground state of SC 3 . 

Moskovits, DiLeila, and Limml25] have carried out resonance Raman 
studies of Ar matrices containing Sc. They report vibrational frequencies 
at 246, 145, and 151 cm -1 which they assign to Sc 3 . Upon warming, 
the two lower frequencies coalesce into one frequency at 150 cm -1 . This 
behavior is interpreted as the removal of a matrix induced assymetry in 
the SC 3 molecule in the originally deposited matrix. As these authors point 
out, the assignment of these lines to SC 3 is based on comparison to other 
systems where the carrier of the spectrum could be established by isotopic 
substitution, but that this is not possible for Sc 3 , due to only one stable 
isotope of the Sc atom. We calculate .symmetric stretch frequencies of 513 
cm -1 for 2 Ag and 204 cm -1 for 2 £^ The large difference in force constants 
between these two molecular states probably results from the steeper inner 
*wall for the state as a result of the larger 4s population. .Assuming 

that is the ground state of Sc 3 , as suggested by a combination of 
the T£SR results and the computations, the large discrepancy between our 
calculated symmetric stretch frequency of 513 cm -1 and the experimental 
value of 246 cm -1 suggests that the experimentally observed frequencies 
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do not arise from a Raman transition which terminates in the 2 A 2 state. 
However, it is possible that the final state is 2 E ; . This would imply that 
the separation between 2 Ai( and 2 E' is very small. The presence of non- 
degenerate bending modes in the final state is consistent with the 2 E ; state if 
it is only weakly Jahn- Teller distorted. The 2 E' state need not be thermally 
populated. However, if both the 2 A^ and 2 E' state were thermally populated 
and the 2 E ; state is only “weakly Jahn-Teller distorted the ESR spectrum 
of the 2 E ; state may be too broad to observe. For Raman transitions the 
2 A5) 'State is connected to ^A^, 2 E', or 2 E ;/ states. Thus, we would not see 
aTtaman transition to the state, but would expectto see transitions 
to the 2 A^ state. These “would be expected to occur near 500 cm -1 , where 
broad absorption is observed[25]. 

The lowest state of Sc^ - is a 3 A^ state which has the configuration: 
3 A^ = 4sa' 1 2 4se' 2 3daj 2 3d^ 2 

This state may be derived from the 2 E' state by removal of one electron from 
the 4se / orbital. For the positive ion the 3d shell is expected to be stabilized 
■with respect to the 4s shell; thus, a configuration with a 3d population near 
-4.0 is more stable for the positive ion, while for the neutral a state with a 
3d population near 3.0 3s more stable. Based on the CASSCF results, the 
ionization potential for removal of one -electron from dhe 4se' orbital of the 
:2 E / state is about 3.3 eVL 

V. Conclusions. 

CASSCF/CCI calculations have T>een presented for Sc^ and Sc^\ The 
^calculations use a ]4s4p3dlf] basis set in -conjunction -with an -effective core 
potential, developed by Hay and Wadt, which replaces the Ne core but 
includes the 3s and 3p core levels along with the valence electrons in the 
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calculations. Molecular tests of the Sc basis set and ECP used here indicate 
results in good agreement with all electron calculations. 

The bonding in the lowest state of SC 3 is found to arise from the three 
ground state(4s 2 3d 1 ) atoms. The 4s bonding involves 4s -*■ 4p promotion 
leading to 4s-4p populations of 4s = 4.55 and 4p = 1.17. The promotion 
energy involved is compensated for by the formation of three metal-metal 
bonds; thus, this bonding mechanism is favorable for the trimer but not 
for the dimer. The 3d electrons pair up to give a doublet ground state in 
agreement 'with the experimental ESR results. The low-lying states of SC 3 
have the configurations: 

2 Aj = 4sai 2 4se ' 4 3dai 1 Sdai , 2 
2 A" — 4SUJ 2 4se ' 4 3daj 2 Sda " 1 
2 E' = 4sai 2 4se ' 3 3da ' 2 3daj > 2 

The 2 E ’ state dissociates to two ground state atoms and one excited state 
atom, while the other two states dissociate to three ground state atoms. 

The ground state is assigned as 2 A^ based on a combination of the 
calculations and ESR experiments, but the 2 E' -state is very low-lying . 
The calculated symmetric stretch frequencies for the 2 A^and 2 E 7 states are 
513 cm -1 and 204 cm -1 , respectively. 

The Jowest state of 5c^~ is a state vrhich las the configuration: 
^A" = 4sa ' 1 2 4se ' 2 3da ^ 2 Sda ^ 2 

This state anay he derived from the 2 E} state by a-emoval of one ^electron 
from the 4se / orbital. Tor the positive ion the 3d shell is expected to be 
-stabilized with respect to the 4s shell; thus^ a state ~wiih a 3d population 
near 4.0 is lower for the positive ion, while for the neutral a state with a 
3d population near 3.0 is lower. 
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Table I.CASSCF Energies Tor Sc 3 (Equilateral Triangle Geometries) 

R(Sc-Sc) 2 1 2 A , 1 2 A" 

4.T5 -137.83351 

5.00 -137.85695 -137.82927 

5.25 -137.86828 -137.8516 -137.85275 

5.50 -137.87083 -137.86332 -137.86572 

5.75 -137.86750 -137.87093 -137.87286 

6.00 -137.86063 -137.87821 
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Table II. CASSCF Energies lor 
8 

R(Sc-Sc) 3 A^ 

5.00 - 137.71671 

5.25 - 137.72811 

5.50 - 137.73042 

5.75 - 137.72665 

6.00 - 137.71903 


Sc^~ (Equilateral 
3e> 


- 137.69792 

- 137.70911 

- 137.71456 

- 137.713 


Triangle Geometries) 

3 E V 


- 137.69379 

- 137.70612 

- 137.71212 


J 


2i ; 


i 



Table III. CCI Energies lor Sc 3 (Equilateral Triangle Geometries) 


R(Sc-Sc) 

5.00 
5.25 
5.50 
5.T5 

6.00 


2 A" 

-137.99087 

-138.01180 

-138.02436 

-138.03413 

-138.01591 


2 E / 


-138.02351 

-138.02899 

-138.02890 

-138.02570 
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Table IV. 8CF and 8CF/CI Energies lor Ca 3 (Equilateral Triangle Geometries) 


R(Ca-Ca) 

SCF 

SCF/CI 

6.00 

-108.5T141 

-108.69048 

6.50 

-108.60954 

-108.72454 

T.00 

-108.62878 

-108.74015 

T.50 

-108.63768 

-108.74544 

8.00 

-108.64148 

-108.74566 

8.50 

-108.64300 

-108.74390 

50.0 

-108.64304 

-108.73066 


23 



617 


Figure 1. Calculated potential curves for the Ca^ molecule for equilateral 
triangle geometries. The calculated binding energy for SCF and SCF/CI 
calculations is plotted as a function of R(Ca-Ca). 

Figure 2. Calculated potential curves for selected states of SC3 and Sc^ - for 
equilateral triangle geometries. As for Fig. 1 the binding energy is plotted as 
a function of R(Sc-Sc). The curves are from CASSCF calculations except for 
the 2 A^ and 2 E / states where contracted Cl curves are also shown. The zero 
of energy for SC3 is taken as the sum of three ground state Sc atoms, while 
for Sc^ the zero of energy is two ground state atoms and one Sc+ atom in 
the 4s 1 3d 1 state. 

Figure 3. Contour plots of selected orbitals for the l e A , 1 state of Ti^~ . The 
orbitals are plotted in the molecular plane except for the 3d7r" orbitals which 
^are plotted 1.0 a 0 above the molecular plane. 

Figure 4. Contour plots of selected orbitals for the l e A / 1 state of Ti^“ . The 
orbitals are plotted in the molecular plane. 
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Calculated Ground State Potential Surface and Excitation Energies 

for the Copper Trimer 
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Abstract 

The results of an SCF/SDCI treatment are presented for selected por- 
tions of the ground state potential energy surface for the C 113 molecule. 
T'or equilateral triangle geometries(D 3 /,) the lowest state is 2 E' arising from 
-4sai 2 4se n . The 2 E ; state exhibits strong Jahn- Teller distortion, leading to 
^Ai (acute angle) and 2 B 2 (obtuse angle) minima in C 2v symmetry. Here 
the 2 B 2 minimum is a true minimum on the surface while the 2 Ai minimum 
is a saddle point connecting adjacent 2 B 2 minima. This strong Jahn- Teller 
distortion is consistent with the observed ground state vibrational levels in 
the gas phase spectroscopic studies of Morse, Hopkins, Langridge-Smith, 
and Smalley. The 2 B 2 minimum is also consistent with the observed ESR 
spectrum of CU 3 in a matrix, which has been interpreted as an obtuse angle 
structure with most of the spin density on the end Cu atoms. A linear 2 £+ 
state is found to be 0.26 eV higher. 

Calculated vertical excitation energies for C 113 indicate that the upper 
=state in the transition studied by Morse et. al. should be ?A' X not 2 E". An 
.^alternativeassignmentibased£on^ :: ?Aj 5 upperrstatelsshown.toiaccount'^or 
^thej>attem of adlowedihot3)ands. 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035 
^Mailing Address: Analatom Inc., 253 Humbolt Court, Sunnyvale, CA 94089 
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I. Introduction 

Transition metals and transition metal(TM) compounds are currently 
of considerable interest, because of their relevance to catalysis and to materials 
science problems such as hydrogen embrittlement and crack propagation in 
metals. For a review of the current state of both theory and experiment for 
TM molecules, the reader is referred to the review article by Weltner and 
Van Zee[l]. Recent progress in the experimental characterization of TM 
molecules has come from both matrix isolation and gas phase spectroscopic 
approaches, while progress in the theoretical treatment of TM molecules has 
occurred as a result of improved methods of treating electron correlation. 
These improvements in part involve more efficient computational methods 
such as improved MCSCF[2],CASSCF[3], and direct CI[4] techniques, but 
also qualitative ideas such as the GVB[5] method have provided a con- 
ceptual framework for understanding electron correlation effects in these 
molecules. 

For the transition metal trimers most of the current experimental in- 
formation comes from matrix ESR studies. Here the assumption is made 
that the electronic state observed in the matrix is the .ground state of 
the molecule. These studies directly give the electronic spin state of the 
molecule, and analysis of the hyperflne structure ;gives some information 
on the electronic distribution of the open shell orbital(s). In some cases 
-limited .geometric information nan also ±>e ^derived, eg. for SC3 all the 
Sc atoms are found to jbe -equivalent, thus implying an -equilateral triangle 
geometry[6]. To date gas phase spectra have been obtained only for the cop- 
per trimer[7]. These spectra permit only a vibrational analysis and thus do 
not give detailed geometric information. However, the vibrational analysis 
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indicates a Jahn- Teller distorted ground state, and provides vibrational fre- 
quencies for an excited state at an excitation energy of about 2.3 eY. These 
new experimental results provide a challenge to theory to understand the 
bonding in the transition metal trimers. 

There have been several previous all electron studies of the CU3 molecule 
[8,9,10]. SCF studies vary as to whether the linear geometry [8] or near equi- 
lateral triangle geometry [9] of CU3 are lowest, while Cl studies [10] using a 
[2sld] contracted basis set, which in effect eliminates most of the 3d shell 
correlation, place the distorted equilateral triangle structures lowest, but are 
mot accurate enough to clearly establish the ground state geometry of this 
molecule. Diatomics in molecules (DIM) studies have also been reported for 

✓ 

CU3 and other elements of group IB of the periodic table[ll]. These studies 
also indicate a near equilateral triangle geometry for CU3. The present 
studies differ from previous studies in that larger basis sets and more exten- 
sive electron correlation are included in order to more accurately define the 
ground state surface. Additionally, calculations of the excited states are 
carried out in order to understand the resonant two photon ionization 
^spectrum oTMorse et. al. [7]. 

In order to understand the bonding in the copper trimer we first con- 
sider some features of the bonding in the copper dimer. Qualitatively, 
the bonding in CU2 has been described as a 4s-4s bond arising out of the 
4s 1 -3d 10 .ground state of the Cnatom[L2]. -Some controversy has arisen as 
to the iole of the 3d electrons in this molecule]13]. 'While the 3d 10 shell of 
Cu would reasonably be expected to remain atomic like since it is a closed 
shell and the 4s orbital is large with respect to the 3d orbital for Cu ( The 
ratio <T4 8 > / <r3d> increases monotonically from 2.03 for Sc to 3.36 


2 



4 


TZ9 


for Cu.[14]); the 3d electrons are not inert, since correlating them shortens 
the bond by 0.19 ao. The resulting bond length is 0.15 ao longer than ex- 
periment; this error is found to result from about equal contributions from 
relativistic effects and higher excitations. Relativistic effects have been 
studied by Martin[15], by Scharf, Brode, and Ahlrichs[16], and by Werner 
and Martin[17], using all-electron treatments and by Laskowski et. al.[18] 
who have added relativistic effects using a relativistic ECP. These authors 
find a relativistic contraction of about 0.08 a 0 . Studies of higher excita- 
tions by Martin[15], by Scharf, Brode, and Ahlrichs[16] and by Werner 
and Martin[17] using CEPA show that the remaining discrepancy is due to 
higher excitations. ( Note that with two 3d 10 shells of "spectator electrons” 
one expects a size consistency problem for CU 2 .) These studies indicate a 
need to correlate the 3d electrons in Cu 3 . Indeed we find a substantial 
bond length contraction in Cu 3 due to correlating the 3d electrons. One 
also expects the same problems for Ci^ as were observed for Cu 2 namely a 
relativistic contraction which would shorten the computed bond lengths and 
size consistency errors. Nonetheless, the non-relativistic SCF plus singles 
doubles Cl model proved to be a useful model for Cu 2 and is expected to 
give the main features of the Cu 3 surface. In a related study Laskowski 
et. al.[18] have carried out calculations for the ground state surface of Cu 3 
which incorporate relativistic effects using a relativistic ECP and use the 
CEPA method [19] to reduce the size consistency problem. 

'Starting with three -4s 1 3d 10 Cu attorns and considering only the 4s 
electrons, leads to a 4sa f 1 2 4se /1 configuration for equilateral triangle geometries 
or a 2 E ; state. The 2 E' state is a maximum on the potential surface 
with Jahn- Teller distortion leading to a splitting into 2 Ai (acute angle) 
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and 2 B 2 (obtuse angle) minima in C 2 V symmetry. Increasing the angle 
relative to the value for the 2 B 2 minimum, leads to a gradual increase 
in the energy, and finally a linear geometry 2 £+ state which is calcu- 
lated to be about 0.26 eV higher. In these calculations we concentrate on 
the C 2 v cuts through the potential energy surface. These calculations 
lead to a calculated well depth (energy lowering for Jahn- Teller distor- 
tion of the 2 E state to the 2 B 2 minimum) and barrier to pseudorota- 
tion (difference in energy between the 2 B 2 minimum and the 2 Ax saddle 
point). 

We also report calculations for the positive ion state ( 4sai 2 ) of Cu^ 
and for low-lying valence and Rydberg states. These calculations permit 
comparison to the experimentally measured ionization potential, and per- 
mit identification of the upper state in the transition reported by Morse, 
Hopkins, Langridge-Smith, and Smalley. [7] 

Section H discusses the basis sets and other technical details of the 
calculations. Section HI discusses the calculated potential surface for the 2 E' 
ground state of CU 3 . Section IV describes the calculations for the excited 
states of Cua- Finally Section V presents the conclusions from this work. 

H. Calculational Details 

The basis set used in these studies is essentially the Wachters’ basis[20] 
for the 2 S state 4s 1 3d 10 of the Cu atom. This basis set is augmented 
by the Wachters’ 4p functions .multiplied bjr 1.5 and the Hay diffuse 3d 
functional]. The inner s and p functions are contracted in a segmented 
fashion using contraction HI as recommended by Wachters, while the 4s and 
4p functions are uncontracted, and the 3d functions are contracted (411) 
leading to an [8s6p3dj basis set. For CU 2 , addition of 4f functions to the 
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basis set decreased R e by about 0.04 ao [12]. This is viewed as a small effect 
compared to the relativistic effects which are being neglected and these 
functions were not included here. The small importance of 4f functions 
for CU 2 is a reasonable result given that the bonding involves mainly the 
4s electrons of Cu, but is in contrast to results for those transition metal 
dimers which involve nd bonding where the 4f functions have a large effect 
[22,23,14]. 

For describing the excited states, it was necessary to add more diffuse 
basis functions, since most of the excited states are Rydberg like. The 
diffuse functions can be added either as atomic centered functions, in which 
case diffuse s and p functions are required on each Cu atom to describe 
s, p, and d Rydberg orbitals, or as molecule centered functions ( at the 
molecular center ), in which case s, p, and d functions are added. The 
latter approach has the advantage of requiring fewer basis functions, and 
also avoids linear dependency problems which can arise if very diffuse 
atomic centered functions are added. Three different basis sets were tried. 
The first, designated as diffuse(l), augmented the valence basis set with s, 
p, and d functions located at the molecular center, which were obtained 
as two term GTO fits(uncontr acted) [24] to hydrogenic like STO 3s, 3p, 
.and 3d functions (f=l/3). The diffuse(2) basis set was constructed like 
diffuse(l) except that a three term GTO fit was used, contracted (21), 
.and one additional .diffuse function uf teach type was added using an -even 
tempered criterion. Finally, diffuse(3) was based on diffuse(2), but omitted 
the most diffuse set of molecule centered functions, and added one set of 
atom centered s and p diffuse functions obtained using an even tempered 
criterion. As discussed in section IV, the diffuse(l) basis set proved adequate 
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for describing the lowest six excited states of Cu3( up to excitation energies 
of about 3.5 eV). 

The integrals were computed in C2,, symmetry using MOLECULE[25]. 
The SCF calculations for the ground state surface were carried out in C 2 v 
symmetry even for equilateral triangle geometries where higher symmetry 
could have been imposed. In the case of the calculations for the excited 
states which were restricted to equilateral triangle geometries, the 

higher symmetry was imposed. The SCF calculations were followed by 
singles and doubles Cl. Here the high lying virtual orbitals corresponding 
to the antibonding combinations of the Is, 2s, and 2p atomic levels were 
removed from the calculation. In general the 4s and 3d electrons were 
included in the calculation (33 electrons for CU3). In each case Davidson’s 
correction[26] was computed and added to the computed Cl energy, in order 
to include the main effects due to quadruple excitations. The SCF and Cl 
calculations were carried out using the SWEDEN system of programs. [27] 

HI. The Ground State Surface 

Table la shows energies for C113 as a function of Cu-Cu distance for 
equilateral triangle geometries. From Table la one sees that Cl leads to a 
bond length contraction for Cu$ as was observed for CU2. Here the Cu-Cu 
bond length is 4.90 ao for SCF and 4.59 ao for Cl, or a contraction of 0.31 
ao compared to 0.21 ao for CU2. Table lb shows similar information for 
-Cu Prom Tables la :and lb one sees that Cu3 .and Cu^ Pave veiy similar 
symmetric stretching mode potential curves. The “vibrational frequencies 
derived from the potential curves are 208 cm -1 for C113 and 214. cm -1 
for Cu^~. This result suggests that the 4se ; orbital of the ground state 
4sa , 1 2 4se' 1 configuration is non- bonding to slightly antibonding. (Note that 
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a difference of 6 cm -1 in vibrational frequency is not significant for the grid 
used here.) Since we identify the upper state in the transition observed by 
Morse et. al. [7] as a Rydberg state (See section IV.), it is appropriate to 
compare the computed Cu^ - frequency to the upper state symmetric stretch 
frequency of 243 cm -1 . The only other reported computed frequency 
for this state is 445 cm -1 from DIM calculations [ 11 ]. There is also an 
experimental value of 354 cm -1 by DiLella et. al. [28], but this experiment 
has been questioned by Morse et. al. [7] and by Weltner and Van Zee [ 1 ]. 

Tables II and HI relate to the Jahn- Teller effect for the 2 E ' state of C 113 . 
These Tables show the variation in energy for the 2 Ai and 2 B 2 symmetries 
as a function of Cu-Cu-Cu angle with two Cu-Cu distances fixed at 4.60 
ao- The same information is also presented graphically in Fig. 1 . From 
Fig. 1 and Tables II and m, one sees that the 2 E' state is a maximum 
on the potential surface with distortion to C 2v symmetry leading to ai 
lowering in the energy for both 2 Ai and 2 B 2 symmetries. The 2 Ai symmetry 
favors acute angle geometries while the 2 B 2 symmetry favors obtuse angle 
geometries. The 2 B 2 surface shows a marked asymmetry leading at larger 
-angles to a linear : 2 JC+ state. Table IV gives energies for linear geometries 
of CU3. This state is found to be about 0.26 eV higher than for the 2 B _2 
minimum. 

Tables V and VI show the energies used in optimizing the geometries 
for .the ^Ai nand ^D^ structures in C 2 V symmetry. The optimal .geometries 
were derived via a two dimensionalparabolic fit to the points given in the 
.tables. From the optimal geometries given in Tables V and VI, one sees 
that the 2 B 2 structure has two short Cu-Cu distances corresponding to two 
Cu-Cu bonds, while the 2 Ai structure has one short Cu-Cu distance and 
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two long Cu-Cu distances. From these geometric parameters one expects 
that the 2 Ai structure may be a saddle point interconnecting adjacent 2 B 2 
minima. This has been confirmed at the SCF level by elongating one of 
the long bonds and simultaneously contracting the other which leads to a 
decrease in the energy. The 2 B 2 minimum on the other hand is expected 
to be a true minimum. 

From the vibronic levels of the ground state, as obtained from analysis 
of the hot bands, Morse, Hopkins, Langridge-Smith, and Smalley[7] con- 
clude that the ground state is strongly Jahn- Teller distorted. Following 
Herzberg [29], for an e' vibration of a 2 E' electronic state, the vibrationless 
level is of 2 E' symmetry while the n=l vibrational level leads to vibronic 
states of 2 E', 2 A', and 2 A^ symmetry. For D 3 h symmetry these levels are 
degenerate, but for a strong Jahn-Telier effect the levels are stabilized with 
the 2 A' and 2 A^ levels still degenerate and below the 2 E' level (See Fig. 2). 
Our interpretation of the spectrum “has The : ?A^ and 2 A^ levels 15 cm -1 
above the "vibrationless level, -while the ^E' level is 100 -cm -1 above the 
vibrationless level. 

Morse et. :al.[7] have computed the 2 E' "vibronic levels for a potential 
with a well depth of 543 cm -1 and a barrier .height of 407 cm -1 . They 
find the 2 A\ level to be only .2.8 cm -1 above the ground vibronic level, 
as compared to an experimental estimate of 15 cm"” 1 . The parameters 
calculated in the present work are a well depth of 555 cm -1 and a barrier 
dieightof 171 'em -1 . It would be interesting to seeif these parameters would 
give better agreement with experiment. 

The 2 B 2 ground state is also consistent with the matrix ESR work of 
Howard, Preston, Sutcliff, and Mile [29], which has been interpreted as an 
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obtuse angle structure with most of the spin density on the end Cu atoms. 

Note that the observed spin density is consistent with a 2 B 2 state, but not 
with a 2 Ai state which would have substantial spin density on the central 
Cu atom. For comparison we show in Table VII. the Mulliken populations 
for the singly occupied b 2 orbital for the 2 B 2 state and the singly occupied 
ai orbital for the 2 Ai state. We should note here that the 2 Ai and 2 B 2 
states are calculated to be close in energy which is consistent with recent 
experimental results for Ag 3 by Kernisant, Thompson, and Lindsay [31], 
where it is shown that matrix effects invert the ordering of these states. 

The calculated binding energies for the Cu 2 and Cu 3 clusters were ob- 
tained from the Cl wavefunctions. Here size consistency problems introduce 
some uncertainty in the computed binding energies. For Cu 2 the super- 
molecule Cl energy (R e = 50. a 0 ) is 0.45 eV above the sum of the Cu atom 
Cl energies. This error results from -the .neglect of quadruple excitations 
which .arise asproducts of double ^excitations from ^each of the Cu.atom 3d 10 
shells. Including Davidson’s correction leads to a supermolecule energy 0.34 
eV below the sum of the atom Cl energies. In these calculations we choose 
to compare the Cl plus Davidson’s correction energy near R e with twice the 
singles doubles Cl energy for the Cu atom. This procedure, while consis- 
tent in the level of excitation, may overestimate the binding energy due to 
the overestimation of the effect of quadruples in the molecular calculation 
by Davidson’s correction. The resultant D e is 1.93 eV compared to an ex- 
perimental value of .2.05 eV[32]. Tor CU 3 a reasonably consistent procedure 
is to compare the Cu 3 2 B 2 Cl plus Davidson’s correction energy to the Cu 2 
Cl plus Davidson’s correction energy plus the Cu atom Cl energy. Note that 
within the limits of the use of Davidson’s correction, we are then compar- 
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ing Cl through quadruples for each case. While hextuple excitations could 
contribute for the trimer, they are expected to be an order of magnitude 
less important than the quadruples. The resultant binding energy is 0.56 
eV for CU 3 relative to CU 2 plus Cu atom, as compared to an experimental 
value of about 1.0 eY [33]. The larger error here may be due to neglect of 
hextuple excitations, and the indication is that single double excitation Cl 
treatments ■will be less accurate as the cluster size increases. 

A question which naturally arises is how accurate are the relative 
energies of the various structural isomers of Cu^. The two major sources 
of error here are size consistency and basis set superposition error. The 
magnitude of the size consistency error can be judged from the effect on 
the R e of Cu 2 . Here size consistency effects are judged to make R e too 
long by about 0.05 a 0 ; from this we conclude that higher excitations are 
imore important as R Cu-Cu decreases. An estimate of the JR. dependence 
of this -effect mear R e unay be jnade from the Cu 2 £>DCI ^potential curve. 
Here compressing the bond by 0.10 ao from R e (the change in the potential 
curve necessary to reduce R e by 0.05 ao) increases the energy by 0.010 eV. 
.Thus it would seem that for small changes in geometry, size consistency 
effects would not substantially distort the potential surface. The other 
potential problem is basis set superposition error. Here we calculate a Cl 
superposition error of 0.06 eV for one atom due to the basis set on an 
adjacent atom, for R = 4.5 ao but relatively slow variation with small 
changes in distance near Tl c . Prom the calculated Cl basis set superposition 
error curve, we estimate an increase in the 2 B 2 ->• 2 Ai separation from 
171 cm -1 to 223 cm -1 and a reduction in the equilateral triangle to 
linear separation from 0.26 eV to 0.14 eY. However, these effects will be 
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compensated for to some extent by the size consistency errors which have 
the opposite phase and are of the same order of magnitude. We conclude 
that a more complete basis set/CI calculation will still find equilateral 
triangle geometries below linear geometries and will find a small 2 Bz -* 
2 Ai separation with 2 B 2 lower. These conclusions have been confirmed by 
the results of Ref. 18 which show a smaller 2 B 2 -*• 2 Ai separation than 
in the present calculations and somewhat shorter Cu-Cu bond lengths (due 
to relativistic effects), but do not change any of the major features of the 
computed ground state surface. The excitation energies will not be effected 
by size consistency effects, since the ground state and excited state have 
similiar geometries and hence the size consistency error is constant. 

Comparing now to other computed results, previous all-electron studies 
[8,9,10] did not include sufficient levels of electron correlation to clearly es- 
tablish the teatures of the_ground -state Jsurface.jUthoughthere was general 
agreement that t he jnear ^equilateral tri angle Tegion of the surface displayed 
a strong Jahn- Teller effect and the 2 Ai and 2 B 2 structures were close in 
energy, the relative energies of the near equilateral triangle geometries and 
linear geometries were not clearly established. DIM [11] studies, which in- 
corporate some of these effects in a semiempiric al fashion through the use 
of accurate potential curves for the diatomics, did lead to energetics in 
reasonable agreement with the present calculations. Here the equilateral 
triangle - linear separation is 0.18 eYin^ood agreement with our best es- 
timate of 0.14 -eV, .and the Jahn-Teller stabilization ^energy is 480 cm -1 
compared to 555cm— 1 in our calculation. However, the computed sym- 
metric stretch frequency is 445 cm -1 , which is much larger than our calcu- 
lated value and, accepting our interpretation of the experiment, also much 
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larger than the experimental value of 243 cm 1 . 
IV. Calculated Vertical Excitation Energies. 
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As discussed in section HI, the ground state of C 113 for equilateral 
triangle geometries is 2 E', arising from a 4sa' 1 2 4se ' 1 configuration. From 
comparison of the calculated symmetric stretching frequencies of CU 3 and 
we see that the 4se' orbital is essentially nonbonding. Given this one 
expects that the low-lying excited states of CU 3 will arise by excitation 
of a single electron from the 4se ; orbital to various virtual orbitals, while 
leaving the 4s 2 configuration intact. This is equivalent to saying that 
the excited states formally have the character of Ttydberg like states, since 
they correspond to a single electron in the field due to Cu^\ The solutions 
to such a problem are expected to be hydrogenic like. The excitation 
energies given in table VHI were obtained from a combination of a valence 
Cl calculation, to determine which low-lying states exist and SCF plus 
singles and doubles Cl calculations to determine the excitation energies. 
Because the excited states of C 113 are Rydberg like, a good approximation 
to the occupied orbitals is obtained from :an improved ^virtual orbital(IVO) 
caIcuIation[34j, whicli in the present case is equivalent to a calculation on 
The positive ion. To determine the mature of The low-lying Rydberg states 
we carried out a valence Cl calculation which consisted of full Cl for the 
three valence electrons among the space of the bound virtual orbitals of the 
positive ion calculation. This calculation leads to the valence Cl values in 
Table VIII. Subsequent “SCF plus singles and doubles Cl was carried out To 
determine the ■excitation energies. 

Figure 3. shows selected valence and Rydberg orbitals for C113 obtained 
from an SCF calculation for Cu^ - . The 4se' orbital which is occupied in 
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the ground state is like a hydrogenic 2p function. This is clearly a valence 
state since the 4se' orbital is largely Cu 4s like. The other component of 
the hydrogenic 2p is an a^ orbital, which turns out to be of mixed valence 
Rydberg character. The lowest excited state of C113 is a state arising 
from this level. This excited state which occurs at 1.25 eV is not dipole 
allowed from the ground state of CU3 without vibronic coupling. 

The remaining five excited states of CU3 which are reported here arise 
from singly occupied orbitals that have the character of n=3 hydrogenic 
levels, i.e. they exhibit Rydberg character. The assignment of these 
levels to n=3 is based on the effective quantum numbers obtained from 
the computed excitation energies and the computed ionization potential. 
These values are consistent with quantum defects of 0.8, 0.2-0.6, and 0.6 
for 3s, 3d, and 3p hydrogenic levels, respectively. It should be noted here 
that the Rydberg n value has no relationship to the principle quantum 
number for the Cu atom. The geometries of these states are expected 
to be very similar to the geometry of Cu^\ From section in. we saw 
that both Cu 3 and Cu^" have D^h minima at nearly the same geometry 
(an equilateral triangle geometry with R Cu-Cu = 4.60 ao). Table VIH 
shows the calculated vertical excitation energies for CU3 at the geometry 
given above. The electronic symmetries here are (with the Rydberg orbital 
character given in parenthesis) 2 Ai(3saj), 2 E'(3de'), 2 E // (3de // ), ^^(Sda'J, 
and 2 E'(3pe'). 

'The excited states discussed so far .arise only from the 4s derived levels 
of the Cu atom. Because the 4s 2 3d e state of Cu atom is only 1.5 eV above 
the 4s 1 3d 10 state it is possible that low-lying states may exist in which one 
Cu atom is in the 4s 2 3d 9 excited state. Such states have been studied by 
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Walch and Goddard [35] for Ni 3 . Here it is found that the central atom 
is 4s 2 3d 8 while the end atoms are 4s 1 3d 0 . The central Ni atom forms sp 
hybrids directed at 180 degrees from each other. These sp hybrids then 
form bonds to the 4s orbitals of the end Ni atoms leading to a linear ground 
state. Similarly for C 113 we find the lowest state of this limit is a 2 A state 
with the 3d hole localized on the central Cu atom (linear geometries). Table 
IX gives computed SCF/CI energies for linear geometries. The computed 
excitation energy is 1.63 eV and the best estimate of the excitation energy is 
1.95 eV (correcting for the error in the 4s 1 3d 10 -*■ 4s 2 3d 8 excitation energy 
of the atom). This excitation energy is for the 2 A state at the optimal 
linear geometry compared to the 2 B 2 component of the ground state at its 
optimal geometry. Given the sharp bending curve for the analagous state 
of Ni 3 we expect that this state will be way above the 2.3 eV excitation 
energy in the transition studied by Morse et. al. [7] for equilateral triangle 
geometries. Additionally, the upper state in that transition has a geometry 
close to that of the ground state which essentially excludes linear states 
from consideration. 

Morse et. al.[7] assign the upper state in the transition which they 
observe at about 2.3 eV as 2 E". This assignment is based on the presence 
of anharmonicity in the upper state vibrational levels in the normal mode 
interpretation, which can be explained by assuming that the upper state 
is 'weakly Jahn-Teller distorted; although, -the upper state levels can also 
he lit to vibrational frequencies of 149 cm -1 -and ;252 cm -1 but with 
significant anharmonicities. The Jahn-Teller distorted upper state model 
requires either a 2 E' or 2 E" upper state. The choice of the 2 E" symmetry is 
consistent with the observation that one vibronic component of the ground 
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state ( 2 A\ or 2 A^ ) is not radiatively coupled to the vibrationless level of 
the upper state. The Jahn-Teller distortion in the upper state is found to 
be very small corresponding to a Jahn-Teller stabilization energy of only 
9 cm -1 . Our initial assumption "was that the upper state •was of 2 E " 
symmetry, and the small Jahn-Teller distortion was consistent with the 
Rydberg nature of the e" orbital in this state. However, the calculated 
excitation energy to this state is 3.32 eV which is about 1.0 eV above the 
observed transition energy. After careful checks on the basis set(see below) 
we conclude that our excitation energy could not be in error by this much. 
Rather, we believe that the 2 A\ state, which we calculate to be at 2.14 eV, 
is the upper state observed at about 2.3 eV. This state also turns out to be 
consistent with the selection rules observed for the hot bands. 

Figure 2 shows the vibronic levels for the ground state of C 113 as 
expected for a strongly Jahn-Teller distorted molecule. The actual locations 
are obtained from analysis of the spectrum of C 113 as measured by Morse 
et. al.[7] (See Figure 4). The main feature of Figure 4 is a vibrational 
progression in the bending mode. Here the 0-+0 band corresponds to band 
-2 and the remaining members of the progression are 5, 7, 9, and 12 (The 
vibrational frequency here is 146. cm -1 .). The strong origin band here is 
consistent with the similar geometries for the two states in the transition. 
The symmetric stretching mode corresponds to 243 cm -1 , and accounts for 
6 and 10. The remaining bands .8 .and 11 Tequire nnharmonicity corrections, 
but can be explained in a normal modes interpretation. A Icey remaining 
feature is the presence of only one hot band corresponding to the 1-+0 
transition (band 1 ), but two hot bands corresponding to the l-*-l transition 
(bands 3 and 4). However, this is readily explained, since for the 2 Ai state 
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the vibrationless level is of 2 A\ symmetry while the n=l level is of 2 E' 
symmetry. The origin of bands 1 and 3 is the 2 E' vibronic level of the ground 
state. Here the separation between 2 and 1 is the same as the separation 
between 3 and 5 and places the 2 E 7 vibronic level of the ground state at 
100. cm" -1 above the lowest vibrational level. Similarly the origin of band 
4 is the 2 Aj or 2 A^ level which is 15cm — 1 above the lowest vibrational level 
based on the separation between 5 and 4. The l-*-0 hot band is absent, 
because transitions from the 2 A! 1 or 2 A^ level to the vibrationless level of 
the upper state are not allowed. The one remaining question is why there 
would be considerable anharmonicity in the upper state levels. From Fig. 
1 one sees that a very large anharmonicity is present in the 2 B 2 component 
of the ground state for the bending mode. This results from the low-lying 
linear configuration. A similar effect is expected for the positive ion, and 
thus anharmonicity in the bending modes of the upper state would not be 
surprising. 

If the ionization potential of CU3 were accurately known experimen- 
tally, it would be easier to calculate the excitation energies of the Rydberg 
states. The best procedure here is to compute energy differences between 
the ion and the Rydberg states, since one expects very little differential 
correlation error in this separation. Unfortunately, for CU3 the I. P. has 
only been bracketed as 4.98 < I. P. < 5.98 by use of diffiferent lasers 
And observing whether ionization occurs J7J. Thus, we had to compute the 
jground state ionization potential by SCF plus singles and doubles Cl for 
the ion and neutral states. The result is 4.90 eV for vertical excitation 
at an equilateral triangle geometry with R Cu-Cu = 4.60 a 0 and 4.97eV 
for abiabatic excitation ( i. e. including the Jahn- Teller stabilization of 
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the lower state). Thus, our calculation indicates that the I. P. of CU3 can 
not be much larger than 4.98 eV. An alternative way to get the ionization 
potential is as the experimental excitation energy to the 2 Ai state plus 
the computed ionization potential of the 2 Ai state. This gives 5.10 eV in 
reasonable agreement with the value above. 

Finally, checks were made of the diffuse basis sets for CU3. Here we 
compared the orbital energy spectrum from a Cu^ - calculation for the three 
basis sets used. The maximum difference in orbital energy was for the Ssa^ 
orbital which was deeper by 0.04 eV and 0.06 eV, respectively, for the 
diffuse(2) and diffuse(3) basis sets, as compared to the diffuse(l) basis set. 
These results indicate that the diffuse portion of the basis set is adequate 
to describe these states to better than 0.1 eV. 

V. Conclusions. 

In this paper we report the results of SCF/SDCI calculations for i)selected 
portions of the ground state surface of CU3 ii) for the ionization potential 
of CU3 and iii)for vertical excitation energies for six low-lying excited states. 

We find the bonding in this molecule to involve the 4s electrons of 
the 4s 1 3d 10 configuration of the Cu atom. The 3d electrons are not inert, 
however, and correlation of the 3d shell leads to significant bond contraction 
effects. 

The ground state of CU3 is found to be a 2 E' state arising from 4sa , 1 2 4se' 1 . 
_A .linear .state is found 4x3 fie -CL26 -cV -higher. The ^E' state ex- 
fiibits strong . Jahn- Teller distortion, leading to ^Ai (acute angle) and 2 B2 
(obtuse angle) minima in C2 W symmetry. Here the 2 B2 minimum is a true 
minimum on the surface while the 2 Ai minimum is a saddle point con- 
necting adjacent 2 B2 minima. These calculations lead to a calculated 
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well depth (energy lowering for Jahn- Teller distortion of the 2 E state to 
the 2 B 2 minimum) of 555 cm -1 and barrier to pseudorotation (difference 
in energy between the 2 B 2 minimum and the 2 Ai saddle point) of 171 
cm - " 1 . From analysis of the hot bands in the experimental spectrum, 
Morse, Hopkins, Langridge-Smith, and Smalley[7] conclude that the 2 A! X 
and 2 A^ levels are 15 cm -1 above the vibrationless level of the ground 
state, while the 2 E ; level is 100 cm -1 above the vibrationless level of 
the ground state. These experimental vibrational levels are consistent 
with a strongly Jahn- Teller distorted e' vibrational level of a 2 E' electronic 
state. 

The 2 B 2 minimum is also consistent with the oberved ESR spectrum 
of C 113 in a matrix, which has been interpreted as an obtuse angle structure 
with most of the radical orbital density on the end Cu atoms. A 2 Ax state, 
on the other hand, would have a large density on the center Cu atom, and 
would not be consistent with the experiment. 

The positive ion of C 113 is 1 A ' 1 arising from 4sa' 2 . The best estimate 
of the ionization potential is 5.10 eY based on the experimental excitation 
energy to the 2 A[ state and the computed I.P. of this state. Experimentally 
this quantity is not well known but has been bracketed to be greater than 
4.98 eV and less than 5.98 eY. The calculated symmetric stretch frequency 
is 208 cm -1 for C 113 and 214 cm -1 for Cu^\ 

The low-lying -excited states of CU 3 -have the character of .Rydberg 
states and require diffuse -basis functions for their description. The upper 
state in the transition observed by Morse et. al.[7] at 2.3 eV has been as- 
signed as 2 E" based on i) significant anharmonicity in the upper state vibra- 
tional levels which can be fit by assuming a weakly Jahn- Teller distorted 
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upper state and ii) observation of selection rules prohibiting the radiative 
coupling of the 2 Aj and 2 A^ vibronic levels of the ground state -with the 
vibrationless level of the upper state. However, based on the calculated ver- 
tical excitation energies we assign the upper state as 2 A' r This assignment 
is also consistent with the observed vibronic selection rules. The presence 
of anharmonicity for at least the bending modes of this state is expected 
due to a low-lying linear configuration which leads to observed large anhar- 
monicity in the 2 B 2 component of the Jahn- Teller distorted ground state. 
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Table la. Energies lor Cu 3 (Equilateral Triangle Geometries) 


R(Cu-Cu) E(SCF) E(CI) 

4.25 -4916.44765 -4917.06800 
4-50 -4Qi6. 46241 -49i7. 07585 
4.75 -4916.46806 -4917.07484 
5.00 A -4916.46856 -4917-06914 
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Table IB. Energies for Cu^~ (Equilateral Triangle Geometries) 


R (Cu-Cu) 

4.25 

4.50 

4.60 

4.75 

5.00 


E(SCF) 

-4916.29195 

-49i6.30T46 

-4916.31072 

-4916.31342 

-4916.31382 


E(CI) 

-4916.88681 

-4916.89562 

-4916.89631 

-4916.89531 

-4916.88999 
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Table II. Energies for CU 3 2 Ai (Bent Geometries. R(Cu-Cu)=4.60 ) 

Cu-Cu-Cu angle E(SCF) E(CI) 

50.0 -4916.44842 -4917.06624 

55.6 -4913.46280 -49i7. 07700 

58.0 -4916.46528 -4917.0TT34 

60.0 -4916.46549 -4917.07622 

70.0 -4916.45808 -4917.06361 

80.0 -4916.44562 
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Table III. Energies for C 113 2 B 2 (Bent Geometries, R(Cu-Cu)=4.60 ) 


Cu-Cu-Cu angle 

50.0 

55.0 

60.0 

65.0 

70.0 

80.0 
90.0 


E(SCF) 

-4916.43485 

-4916.4^595 

-4916.46538 

-4916.46965 

-4916.4715^ 

-4916.47235 

-4916.47198 


E(CI) 


-49i7. 06988 
-4917.07634 
-4917.07802 
-4917.07781 
-4917.07587 
-4917.07406 
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Table IV. Energies lor Cu 3 (Linear Geometries) 


R(Cu-Cu) 

4.30 

4.60 

4.90 


E(SCF) E(CI) 

-4916.46000 -4917.06744 

-4916.46819 -49i7- 06943 

-4916.46895 -4917.06443 
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Table V. Geometry Optimization lor CU 3 2 Ai 


Ri 

r 2 

4.60 

4.60 

4.75 

4.60 

4.90 

4.60 

4.60 

4-248 

4.60 

4.460 

4.75 

4.30 


E(8CF) 

-49i6. 46599 
-49i6. 46869 
-4916.47076 
-4916.46280 
-4916.46528 
-4916.46682 


E(CI) 

-4917.07622 

-4^17.076di 

-4917.07627 

-4917.07700 

-4917.07734 

-4917.07790 

-4917.07809“ 
six points. 


4.753 4.383 

0 Predicted minimum from tbe above 
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Table 

VI. Geometry Optimization for 

Cu 3 2 B 2 

Ri 

r 2 

E(SCF) 

E(CI) 

4.60 

4.60 

-4916.46538 

-4917.07634 

4.60 

4.943 

-4916.46965 

-4917.07802 

4.60 

5.277 

-4916.47153 

-4917.07781 

4.50 

4.943 

-4916.46855 

-4917.07876 

4.40 

4.943 

-4016.46647 

-4917.07858 

4.4T 

5.03 

-4916.46891 

-4917.07886 

4.4T4 5.035 -4917.07887 

a Predicted minimum from the above six points. 
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Table VII. Radical Orbital Mulliken Populations 



2 Ai 


CM 

TO 

CM 



Center 

End 

Center 

End 

s 

.63237 

.22106 

.00000 

.83775 

p 

.00272 

. 14392 

.11584 

.04641 

d 

.00103 

.00096 

.00092 

.00093 
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Table VIX1. Calculated Vertical Excitation Energies for Cu 3 . 


State Sym. 

Valence Cl 

SDCI 

2 E'(3pe') 

2.92 


2 k[ Od^) ■ 

2.88 


2 E" (3de") 

2.80 

3.32 

2 E'(3de') 

2.14 


2 A' Osa* ) 

1.79 

2.14 

2 A^(2pa^) 1 

1.06 

1.25 

2 E , (2pe / ) 

0.00 

0.00 
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Table IX. Energies lor C 113 2 A State (Linear Geometries) 


R(Cu-Cu) 

4.30 

4.60 

4.90 


E(SCF) 
-4916. 44501 
-49i6. 45574 
-4916.45630 


E(CI) 

-491T. 01655 
-49iT. 01898 
-4917.01226 
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Figure 1. Calculated bending curves for two C 2 v cuts through the ground 
state potential surface of CU3. Two sides are fixed at R Cu-Cu = 4.60 ao 
and the Cu-Cu-Cu angle is varied. 

Figure 2. The vibronic levels for an e' vibration of a 2 E' molecule with D^h 
symmetry for i)a symmetric state and ii) a strongly Jahn- Teller distorted 
state (from Ref. 29). The locations of the vibronic levels for the Jahn-Teller 
distorted state correspond to the observed ground state vibronic levels for 
CU3 derived from the spectroscopic studies of Morse, Hopkins, Langridge- 
Smith, and Smalley (Ref. 7). 

Figure 3. Contour plots for selected valence and Rydberg orbitals of CU3 . 
Solid lines indicate positive contours, dotted lines indicate negative contours, 
and nodal surfaces are indicated by dashed lines. 

Figure 4. The experimental spectrum of CU3 observed by Morse, Hopkins, 
Langridge- Smith, and SmaIIey(Ref. 7). 
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